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1.0  INTRODUCTION 

In  view  of  the  increasing  role  of  lasers  for  current  and  potential  military  applications, 
information  on  the  propagation  characteristics  of  a laser  beam  through  various  environments 
is  essential.  Although  a wide  variety  of  phenomenon  affects  the  propagation  of  laser  beams, 
only  the  combined  influence  of  absorption  and  scattering  have  been  measured  in  the  current 
experiments.  The  power  levels  were  sufficiently  low  that  no  thermal  blooming  or  plasma 
generation  was  encountered.  The  test  environment  was  limited  to  snow  crystals  falling 
through  air,  and  measurements  were  made  at  wavelengths  of  both  10.6  and  0.6328  n 
for  direct  comparisons. 

Experimental  measurements  of  laser  propagation  through  hard-to-define  mediums  such 
as  snow  are  especially  valuable,  in  view  of  the  analytic  difficulties  encountered  in  modeling 
a snow  environment.  The  experimental  measurements  can  provide  check  points  for 
analytical  models  and  improve  the  estimation  of  the  field  performance  of  laser  devices. 

2.0  EXPERIMENTAL  DATA 

2.1  EXPERIMENTAL  FACILITY 

The  experimental  arrangement  is  shown  in  Fig.  1.  The  equipment  was  installed  in 
a 60-ft-long,  8-ft-diam,  steel  tank  which  can  be  closed  and  evacuated  to  any  desired  pressure 
between  an  atmosphere  and  a few  mm  Hg.  For  the  present  tests,  the  pressure  was 
atmospheric. 

Both  a CO2  laser  HeNe  laser  were  used  concurrently  in  the  experiment,  so  that 
simultaneous  attenuation  measurements  could  be  made.  The  lasers  were  installed  at  one 
end  of  the  tank,  with  the  CO2  beam  slightly  above  the  HeNe  laser  beam.  At  the  other 
end  of  the  tank,  the  detectors  were  installed  with  the  CO2  detector  below  the  HeNe 
detector.  This  arrangement  made  the  two  laser  beams  essentially  coincident  through  the 
intervening  snowfield,  so  that  attenuation  measurements  for  the  two  beams  are  directly 
comparable. 

The  CO2  laser  was  a continuous  three-watt  device.  The  initial  beam  was  expanded 
to  a 1.7-in.  diameter,  using  a germanium  lens,  and  then  focused  onto  the  detector  at 
the  far  end  of  the  tank.  This  focusing  was  accomplished  with  a 6-in.-diam,  48-in.  focal 
length  spherical  concave  mirror.  The  beam  size  at  the  detector  was  approximately  1/4 
in.  in  diameter.  The  measured  power  delivered  to  the  detector  was  0.7  S watt. 
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The  first  part  of  the  test  was  conducted  using  a detector  manufactured  by  Molectron 
Corporation  to  measure  the  CO2  beam  power.  This  detector  required  the  use  of  a chopper 
to  modulate  the  laser  beam,  and  additional  electronics  to  electronically  remove  the 
chopping  perturbations  from  the  output  signal. 

During  the  latter  part  of  the  test,  a detector  manufactured  by  the  Santa  Barbara 
Research  Center  was  used  which  did  not  require  the  use  of  a chopper  and  the  associated 
electronics.  This  detector  required  liquid-nitrogen  cooling.  The  lasers  and  detectors  used 
are  listed  in  Table  1. 

The  HeNe  laser  used  was  a 15-mw  continuous  device.  The  beam  was  expanded  to 
1.5  in.  in  diameter  and  then  focused  onto  the  detector  at  the  far  end  of  the  tank.  The 
expanding  and  focusing  were  accomplished  using  a commercial  beam  expander  which 
attached  directly  to  the  laser. 

At  the  snowfield,  both  beams  were  approximately  1 in.  diameter  and  coincident. 

The  snowfield  generator  used  is  28  in.  long,  12  in.  high,  and  is  mounted  24  in. 
directly  above  the  laser  beams,  aligned  along  the  beams.  This  produces  a snowfall  which 
is  about  3 in.  wide  and  28  in.  long,  which  falls  through  the  coincident  laser  beams.  The 
snowfield  generator  and  operation  are  described  in  Ref.  1. 

Figure  2 is  a sketch  relating  the  snowfield,  laser  beam,  and  field  camera  viewing 
direction.  A still  camera  is  positioned  to  view  the  snowfield  from  a direction  perpendicular 
to  the  laser  beam  as  indicated  in  the  sketch.  Illumination  is  provided  by  strobe  li^ts 
to  give  a stop-action  picture  of  the  snowfield  for  posttest  analysis. 

2.2  TEST  PROCEDURE 

Starting  with  the  lasers  aligned  and  operating,  the  snow  generator  plate  is  cooled 
to  operating  temperature  (-90°F)  and  maintained  at  this  temperature.  The  plates  are  then 
cleared  by  knocking  off  the  accumulated  frost. 

The  snowfall  density  is  controlled  by  varying  the  snow  growth  time.  For  the  present 
tests,  it  was  found  that  measurable  laser  attenuation  could  be  obtained  using  snow  growth 
times  in  excess  of  two  minutes,  with  a ten-minute  growth  time  completely  attenuating 
the  signal  to  the  detectors. 

Once  the  desired  growth  time  has  lapsed,  the  snow  generator  is  impacted  with  a 
solenoid  plunger  to  initiate  the  snowfall.  After  a time  delay  of  0.5  sec,  transient  recorders 
are  activated  to  record  the  detector  outputs.  Typically,  the  snowfall  begins  to  arrive  at 
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the  laser  beam  at  about  0.7  sec,  with  the  last  of  the  snowfall  passing  through  the  beam 
at  about  1.1  sec.  At  0.9  sec,  strobe  lights  are  flashed  to  illuminate  the  snowfield  for 
the  field  camera  picture. 

After  each  test  the  contents  of  the  transient  recorders  are  displayed  on  oscilloscope 
screens  and  recorded  photographically.  The  laser  beams  are  then  blocked  to  establish  the 
complete  attenuation  levels,  and  these  levels  are  also  displayed  on  the  respective  oscilloscope 
screens  and  recorded  on  the  same  photographs.  A typical  test  result  is  shown  in  Fig. 
3.  Note  that  the  position  of  the  centerline  of  the  laser  beam  at  0.9  sec  is  indicated  in 
the  snowfield  picture.  The  snowflakes  along  the  centerline  produced  the  laser  beam 
attenuation  recorded  at  0.9  sec.  Prior  attenuation  resulted  from  the  interaction  with  the 
snowflakes  which  had  previously  crossed  the  indicated  laser  beam  centerline,  while 
subsequent  attenuation  is  caused  by  the  snowflakes  recorded  above  the  laser  beam 
centerline  in  the  snowfield  picture.  Similar  test  results  for  all  test  data  are  presented  in 
Appendix  A. 

The  test  log  (Table  2)  lists  the  pertinent  parameters  for  the  tests. 

2.3  SNOWFIELD  CHARACTERIZATION 

The  posttest  characterization  of  the  snowfield  is  accomplished  from  analysis  of  the 
single  field  camera  picture  obtained  on  each  run.  The  location  of  the  laser  beam  centerline 
is  established  in  the  snowfleld  picture.  The  field  is  subdivided  into  3-  by  3-in.  areas,  with 
the  first  strip  centered  along  the  laser  beam  centerline  as  indicated  in  Fig.  2.  Subsequent 
strips  of  squares  are  located  at  intervals  of  3 in.  above  and  below  the  centerline. 

Each  3-  by  3-in.  area  is  examined  using  a Quantimet®  720  image  analyzing  computer 
which  counts  the  number  of  particles  in  a preselected  size  range  and  gives  the  total  area 
of  particles  in  that  size  range  contained  with  the  3-  by  3-in.  area.  Twenty-two  different 
size  ranges  are  considered.  This  information  is  recorded  on  paper  tape  for  eventual  computer 
processing. 

The  3-in.  size  represents  the  lateral  depth  of  the  snowfleld.  Therefore,  a side  view 
of  a 3-  by  3-in.  area  contains  all  the  flakes  included  in  a 3-  by  3-  by  3-in.  cell.  The 
population  information  for  each  3-in.  cell  in  the  snowfleld  is  recorded  and  available  for 
processing.  A sample  of  the  available  snowflake  size  distribution  data  is  given  in  Table 
3 for  a particular  cell  on  a particular  run.  These  raw  data  are  available  for  each  ceU  in 
each  snowfleld. 

These  snowfield  data  are  computer  processed,  starting  with  the  uprange  cell  through 
which  the  laser  beams  enter  the  snowfall.  The  total  area  of  snowflakes  within  each  size 
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range  is  divided  by  the  number  of  flakes  in  that  size  range  to  find  an  average  snowflake 
area.  This  area  is  assumed  to  be  circular,  and  the  resulting  diameter  is  then  found.  This 
diameter  is  used  for  subsequent  calculations. 

Snowflake  number  density  is  assumed  to  be  uniformly  distributed  within  a cell,  but 
may  differ  from  cell  to  cell.  The  assumption  of  uniform  distribution  within  a cell  is 
necessary  to  infer  the  number  of  flakes  that  populate  the  l-in.-diam  cylindrical  volume 
of  the  laser  beams  within  that  cell  at  a particular  time. 

The  snowflakes  along  the  laser  beam  centerline  in  the  field  camera  picture  produce 
the  measured  attenuation  at  0.9  sec.  The  cell  contents  are  also  analyzed  for  cells  located 
in  intervals  of  3 in.  above  and  below  the  centerline.  It  is  assumed  that  the  snowfall  velocity 
is  such  that  0.08  sec  give  a fall  of  3 in.  Therefore,  sampling  the  cells  3 in.  above  the 
centerline  gives  the  snowfield  distribution  at  0.92  sec,  while  the  cells  centered  3 in.  below 
the  centerline  give  the  distribution  at  0.82  sec.  This  correspondence  between  the  snowflake 
distribution  at  some  elevation  at  time  equal  0.9  sec  and  the  distribution  on  the  laser 
beam  centerline  at  another  time  assumes  a uniform  falling  velocity  for  all  size  snowflakes. 
Since  large  flakes  fall  faster  than  smaller  flakes,  this  correspondence  becomes  increasingly 
inexact  with  increasing  distance  from  the  centerline  at  0.9  sec. 

2.4  ATTENUATION  DATA  REDUCTION 

The  laser  beam  attenuation  is  measured  directly  from  the  output  of  the  detectors. 
The  no-snow  case  is  taken  to  give  complete  signal  transmission,  while  the  complete 
attenuation  detector  signal  level  is  determined  by  completely  blocking  the  laser  beam. 
The  detector  output  is  a linear  function  of  the  laser  beam  transmission. 

The  individual  effect  of  large  flakes  can  be  seen  in  the  early  portions  of  the  data 
traces,  for  both  laser  wavelengths.  Also,  the  CO2  detector  gave  a noisy  signal,  so  that 
some  smoothing  is  required  to  extract  the  attenuation  data.  The  data  prints  shown  for 
comparison  with  the  calculated  attenuation  values  were  read  from  the  traces  at  selected 
time  intervals.  The  attenuation  data  are  tabulated  in  Table  4,  which  gives  the  ratio  of 
the  attenuated  signal  to  the  no-snow  signal  levels,  as  a function  of  time. 

3.0  ANALYSIS 


3.1  INTRODUCTION 

This  analysis  is  intended  to  aid  in  the  interpretation  of  the  present  laboratory 
experiments.  The  analysis  is,  therefore,  specifically  directed  toward  this  limited  objective. 
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The  analysis  has  been  divided  into  three  parts,  the  first  of  which  consists  of  evaluating 
the  effect  of  a single  scattering  particle.  The  second  part  models  the  effect  of  multiple 
scattering  particles,  and  in  the  third  part,  the  results  are  applied  to  the  experimental 
snowfields  to  give  laser  beam  attenuation  predictions  for  the  present  tests. 

3.2  ATTENUATION  BY  SINGLE  PARTICLE 

The  experimental  arrangement  assumed  has  been  shown  in  Fig.  1 , and  the  pertinent 
dimensions  are  given  in  the  sketch  (Fig.  4).  The  laser  beam  is  assumed  to  be  of  uniform 
intensity,  collimated,  and  be  at  least  as  large  as  the  detector  diameter.  The  detector  diameter 
is  assumed  to  be  larger  with  respect  to  the  scattering  particle.  The  particle,  in  turn,  is 
assumed  to  be  much  larger  than  the  laser  wavelength.  Finally,  the  distance  from  the  particle 
to  the  detector  (fi)  is  assumed  to  be  large  with  respect  to  any  other  dimensions. 

The  particle  is  assumed  to  be  a sphere,  and  under  the  present  set  of  assumptions, 
Fraunhofer  differentiation  theory  can  be  employed  (Ref.  2).  Figure  5 depicts  the  radiation 
intensity  pattern,  which  is  anticipated  on  the  detector.  First,  all  the  incident  radiation 
in  an  area  of  twice  the  geometrical  cross  section  of  the  particle  is  scattered.  Therefore, 
only  scattered  radiation  appears  within  this  area  on  the  detector  surface.  Outside  this 
area,  both  the  incident  and  scattered  radiation  are  present  and  are  assumed  to  be  additive. 


The  power  to  the  detector  is  given  by 
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Then, 


= -2rr 


y\ 


(2) 


Ac  is  the  effective  area  of  extinction.  The  power  to  the  detector  is  then  given  by 

A 


This  is  equivalent  to  assuming  the  initial  intensity  (lo)  acts  on  the  fraction  of  detector  area 
(1  - Ac/Al),  while  no  radiation  is  received  by  the  fraction  of  the  detector  area  given  by 
Ac/Al. 

The  equation  expressing  Fraunhofer  diffraction  intensity  distribution  (Ref.  3)  is 


= _L  v4  r 

k2  ^ L X sin  0 J 


(3) 


where 


X = ;rDp/A 
k2  = (277  J A)  2 

and  Ji  is  the  first-order  Bessel  function. 

The  ratio  of  the  detector  signal  expected  with  a single  scattering  particle  to  that 
expected  without  scattering  is 


_p 

p 

O 


(4) 


where 


O 


(5) 


Equation  (5)  has  been  evaluated  for  both  0.6328  and  10.6  n radiation,  using  fi  = 
28  ft  and  = 0.5  in.  The  results  shown  in  Fig.  6 compare  these  calculations  and  the 

results  with  no  scattering  in  the  direction  of  the  detector.  For  the  shorter  wavelength 
radiation,  scattering  in  the  direction  of  the  detector  makes  an  appreciable  contribution 
to  Ae. 
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The  results  in  Fig.  6 can  be  approximated  by 

A = exp  j -19.124707  + 1.8015599  in  D„  + 0.00539783  (JnDJ^l  (6) 

0.6328  P P ’ 

A = exp  1-21.2492  + 2.491378  in  - 0.04095798  (inDJ^l  (7) 

^106  ' ° P* 

with  Ae  in  in.2  and  Dp  in  microns,  and  where  fin  Dp  is  the  natural  log  of  the  particle 
diameter. 


Figure  7 shows  the  signal  attenuation  due  to  a single  scattering  particle.  The 
contribution  of  the  forward  scattering  component  is  frequency  dependent,  and  the  results 
for  0.6328  and  10.6  fi  are  compared.  The  attenuation  is  shown  to  increase  with  both 
particle  size  and  wavelength.  This  comparison  is  specifically  for  fi  = 28  ft,  D^  =0.5  in., 
and  Dl  = 2 in. 

3.3  ATTENUATION  BY  MULTIPLE  PARTICLES 

The  experimental  configuration  is  such  that  the  laser  beam  is  attenuated  by  a cloud 
of  particles.  This  configuration  has  been  depicted  in  Fig.  2.  The  field  has  a limited  width 
(W)  and  in  practice  is  viewed  from  a direction  perpendicular  to  the  path  of  the  laser 
beam  by  a field  camera.  The  field  camera  is  used  for  posttest  evaluation  of  the  attenuation 
environment. 


The  field  is  subdivided  into  cubes  of  dimension  W for  analysis,  as  shown  in  Fig. 
8.  The  side  view  gives  the  number  of  particles  in  the  various  size  ranges  contained  within 
the  cell.  The  number  of  particles  in  the  size  range  in  the  laser  view  is 


^L.i  - '^SV.i 


4W^ 


(8) 


The  particles  within  a particular  cell  are  assumed  not  to  interact.  Therefore,  the 
scattering  from  various  particles  in  a common  cell  is  taken  to  be  additive.  The  radiation 
which  is  transmitted  through  a cell  is 


P = I. 


M 

' i=l  ' 


(9) 


where  Agj  is  the  extinction  area  of  a sin^e  F**  sized  particle;  Nlj  is  the  number  of 
particles  in  the  it**  size  range  as  seen  by  the  laser  beam;  and  M is  the  number  of  size 
increments  considered  within  the  cell. 


The  extinction  area  is  given  by  Eq.  (6)  or  Eq.  (7),  and  the  Nsv,i  values  are  found 
from  analysis  of  the  field  camera  pictures.  The  radiation  incident  on  the  first  cell  is  assumed 


11 


AEDC-TR-77-99 


to  be  li,  acting  on  area  Al.  After  attenuation  by  the  first  cell,  the  radiation  incident  on  the 
second  cell  is  assumed  to  be  Iq  incident  on  Al  (1  - Aj/Al)  with  no  radiation  incident  on 
area  Ai . Ai  is  defined  as 


A. 

J 


M 


N,  - A 
L,i  e. 


cell  I 


(10) 


with  j = 1. 


The  radiation  incident  on  subsequent  cells  is  reduced  by  a factor  of  (1  - Aj/Alj^) 
for  each  cell  which  is  traversed.  The  radiation  incident  on  the  Kt**  cell  is,  therefore,  lo  on 
an  area  Ai^,k  with  no  radiation  on  (Al  - Ai^,k)-  is  given  by 


K-  1 

^1  .K  = ^ 

0 ,=  1 


(11) 


This  model  is  depicted  in  Fig.  9.  In  summary,  no  interactions  are  considered  between 
particles  in  a common  cell,  but  interactions  are  allowed  between  cells.  By  using  this  model, 
the  signal  after  attenuation  by  K cells  relative  to  the  no  attenuation  signal  is 


3.4  LASER  ATTENUATION  PREDICTIONS 


(12) 


The  well  established  analysis  for  scattering  by  a sphere  has  been  applied  in  an 
approximate  manner  to  the  case  of  multiple  scattering  particles.  The  scattering  performance 
of  real  particles  with  non-spherical  shapes  is  assumed  to  be  approximated  by  the 
performance  of  ideal  spherical  shapes. 

The  spherical  shapes  are  described  in  terms  of  the  diameter.  The  scattering 
performance  of  other  shapes  may  be  adequately  correlated  by  some  effective  diameter 
which  is  related  to  the  particle  dimension.  In  the  absence  of  other  data,  the  mean  particle 
diameter  is  taken  to  be  the  effective  spherical  diameter.  The  mean  particle  diameter  is 
the  square  root  of  the  particle  area  in  the  field  camera  picture. 


The  calculation  procedure  has  been  applied  to  each  snowfield  in  the  present  tests, 
for  various  times  during  the  snowfall  for  comparison  with  the  measured  attenuation. 

The  individual  contribution  of  each  size  particle  in  the  cell  is  calculated  according 
to  Eq.  (6)  or  (7),  weighted  by  the  number  of  particles  of  that  size  range,  using  Eqs. 
(8)  and  (10),  and  the  final  result  is  obtained  from  Eq.  (12).  A sample  comparison  is 
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shown  in  Fig.  10  for  the  snowflake  population  encountered  on  runs  9 and  10.  Comparisons 
are  shown  in  Appendix  B for  all  other  tests. 

Agreement  between  predicted  and  measured  attenuation  would  imply  that  the  random 
shaped  snowflakes  scatter  radiation  with  the  same  efficiency  as  spheres  with  the  same 
diameter  based  on  matching  areas.  The  experimental  results  indicate  that  the  scattering 
efficiency  of  a snowflake  may  be  slightly  larger  than  for  a sphere  vdth  the  same 
cross-sectional  area. 

Since  the  intervening  snowflake  population  is  best  defined  at  0.9  sec,  the  comparisons 
between  theory  appropriate  to  spherical  particles  and  the  experimental  attenuation  have 
been  made  at  0.9  sec.  The  results  are  shown  as  Figs.  11  and  12,  which  are,  respectively, 
appropriate  to  the  0.6328-  and  10.6-/!  wavelengths.  The  agreement  between  the  calculated 
and  measured  attenuation  is  probably  within  the  experimental  data  scatter.  The  best 
quadratic  fit  is  also  shown  in  the  figures,  but  it  is  not  apparent  that  this  representation 
of  the  data  is  more  justified  than  taking  the  calculated  attenuation  to  represent  the  data. 

3.5  LASER  ATTENUATION  PREDICTION  FOR  NATURAL  SNOWFALL 

The  calculation  model  developed  has  been  supported  by  the  experimental  data.  It 
would  be  of  interest  to  calculate  the  attenuation  of  a laser  beam  through  a natural 
enviromnent. 

If  the  ranges  are  long  so  that  the  angle  subtended  by  the  target  is  small,  some 
simplifications  can  be  made  in  the  calculations.  For  instance,  if  the  range  is  2,000  ft 
and  the  laser  beam  at  the  target  is  1 in.  in  diameter,  then  the  an^e  subtended  by  the 
target  is  small,  and 


^max  = " 1-19  X 10"^  deg  (13) 

This  small  angle  allows  the  intensity  function  (Is/Io)  to  be  approximated  by  the  value 
^max  0)  i-®-» 


I 

B 

I 


0 


lim  X*  rjj(X  sin 
0-*o  . X sin0  _ 


Jil 

4k2 


(14) 


The  extinction  area  for  a single  particle  is  given  in  Eq.  (2),  which  using  Eq.  (14)  becomes 


A 


e 


dd 


(15) 
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where 


2 nr  X* 

4k  2 


and 


do  = 1 


9 


Ian  0, 


The  extinction  area  for  a single  particle  is,  therefore,  given  by 


rrD^  r 

1 - 1 

/ "■1-^  V 
— ) 

tan^  0 1 

2 L 

8 

V A / 

m ax  1 

(16) 


(17) 


(18) 


This  expression  is  limited  to  small  values  of  0m  ax  and  replaces  Eqs.  (2),  (3),  and  (5) 
with  this  restriction. 

The  next  step  is  to  extend  this  result  for  scattering  from  a single  particle  to  the 
case  of  multiple  particles.  This  is  done  using  Eq.  (9),  which  assumes  no  interaction  between 
particles  in  a common  cell. 


P = 1. 


The  total  extinction  area  of  a cell  is  given  by  Eq.  (10), 


M 

A.  = V'  N,  ■ A 
i=  1 


cell  j 


(19) 


(20) 


where  the  summation  is  performed  over  the  number  of  snowflakes  in  the  cell  in  the 
i**’  size  range  within  the  laser  beam.  Interaction  between  cells  is  given  by  Eq.  (12), 


(21) 
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Substitution  of  Eq.  (20)  into  Eq.  (21)  gives: 


p 


(22) 


where  Al  is  the  cross-sectional  area  of  the  laser  beam  in  the  j*'’  cell  and  K is  the  number 
of  cells  traversed  by  the  laser  beam. 


The  number  of  flakes  in  a given  size  range  within  the  laser  beam  in  a particular 
cell  (NL,i)  is  related  to  the  unit  number  density  of  flakes  in  space  (N|)  and  the  cell 
dimensions  by 


N 


L.,i 


N: 


A. 


AD 


p,i 


(23) 


where  Lc  is  the  extent  of  the  cell  in  the  direction  of  the  laser  beam  and  ADp.  is  the 
particle  size  range  increment  for  the  i***  size  range.  The  upper  product  limit  (K)  in  Eq. 
(22)  is  determined  by  the  total  extent  of  the  snowfleld, 


K 


= L.,/L, 


Equation  (23)  becomes,  with  these  substitutions, 


(24) 


_P. 

P 

0 


n 

j=i 


(25) 


If  all  cells  are  assumed  to  have  a similar  snowflake  size  distribution,  then  Eq.  (25) 
can  be  replaced  by 


_p. 

p 

0 


[ 


M 

S ad 

i = 1 


P.i 


N: 


L„/L 

I C 


(26) 


The  appropriate  cell  length  (Lc)  is  not  known,  but  a 3-in.  length  has  been  used  for 
the  analysis  of  the  present  experimental  data.  If  the  cell  length  approaches  zero,  Eq.  (26) 
becomes 


lim 

L -»0 
c 


Lt/L, 


exp 


AD 


(27) 
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The  final  equation  for  the  attenuation  of  a laser  beam  due  to  a natural  snowfall  is,  therefore, 

M 

= exp  ADp  .N.  ,A^ (28) 

where  P/Pq  is  the  ratio  of  the  laser  power  exiting  the  snowfall  to  the  laser  power  incident 
on  the  snowfall,  S is  the  distance  into  the  snowfall,  ADp^i  is  .the  particle  size  range 
increment,  N,  is  the  total  number  of  particles  in  the  P*'  size  range  within  a unit  volume 
per  unit  size  range,  and  Aej  is  the  extinction  area  due  to  the  presence  of  a single  particle 
in  the  size  range.  The  analysis  is  restricted  to  large  distances  from  the  snowfield  to 
target  resulting  in  small  angles  (dm  ax)  subtended  by  the  target.  The  quantity  (A^j)  is 
obtained  from  Eq.  (18). 

Information  on  the  number  density  of  snowflakes  in  a natural  snowfall  and  an 
equation  which  expresses  the  number  density  as  a function  of  liquid  water  content  and 
particle  dimension  was  obtained  from  AFSC/SAMSO: 

N;  = 0.323  (a/D„)*  B exp  (.a  D/D^)  (29) 

where  Nj  is  the  number  of  particles  per  unit  volume  per  unit  size  range  increment  of 
melted  diameter,  and  Do  is  an  "equivalent  melted  value"  which  is  a function  of  the 
liquid  water  content  (m)  of  the  cloud. 

= 0.83  m°-286  (30) 

for  small  snow,  and 

= 1.17m®-226  (31) 

for  large  snow. 

The  term  Dq  is  in  millimeters  with  m in  g/m^,  p is  the  density  of  water,  taken 
to  be  1 g(cm? , and  a is-  a constant  3.67.  The  equivalent  melted  diameter  (Dj)  is  related 
to  the  immelted  dimension  by 

D.  = 0.4  D ° 

I P 

D.  = 0.37  D 
* P 

for  small  snowflakes,  and 


D < 0.5  MM 

(32) 

Dp  .>  0.5  MM 
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Dj  = 0.4  Dp®''®2 
D.  = 0.4  D 

‘ P 

for  large  flakes. 

Figure  13  shows  the  numerical  evaluation  of  Eq.  (28)  for  the  case  of  dm  ax  = 

X 10'^  deg.  The  correlation  of  flake  type  (large  or  small)  versus  liquid  water  content 
(m)  has  been  taken  from  a representative  storm  as  given  in  Ref.  4.  From  Fig.  13,  the 
attenuation  due  to  1,000  ft  of  a cloud  of  large  snowflakes  with  a liquid  water  content 
of  0.2  g/m^  is 


Dp  < IMM 
D > IMM 

p 


(33) 


_L  = exp  (-1000  x0.00177)  = 0.17 
P 

O 

Thus,  the  attenuation  of  a CO2  laser  beam  by  a snowfall  appears  to  be  a signiflcant  problem. 

From  the  present  study,  it  can  be  concluded  that  the  attenuation  of  a low  power 
density  laser  beam  by  a snowfall  can  be  analyzed  using  geometric  optics.  For  laser  beams 
with  power  densities  sufflciently  large  to  effect  the  snowflakes,  the  laser  attenuation  may 
be  appreciably  different. 


4.0  RESULTS  AND  CONCLUSIONS 

The  simultaneous  attenuations  of  a CO2  and  a visible  HeNe  laser  beam  by  a common 
snowfleld  have  been  measured  in  a laboratory  experiment.  The  snowflelds  encountered 
have  been  characterized  by  counting  the  number  of  flakes  in  various  size  ranges. 

In  conjunction  with  the  experiment,  a calculation  model  has  been  developed  to  predict 
the  attenuation  due  to  individual  flakes.  This  model  assumes  that  the  flake  behaves  as 
a sphere  with  the  same  cross-sectional  area.  The  particle  intercepts  the  radiation  from 
an  area  of  twice  the  particle  geometric  area,  but  some  fraction  of  this  radiation  is  scattered 
in  the  direction  of  the  target. 

The  experimental  laboratory  results  support  the  calculation  model,  which  predicts 
a wavelength  dependence  on  the  quantity  of  radiation  scattered  toward  the  target,  but 
no  such  dependence  on  the  radiation  intercepted  by  the  particle.  Therefore,  as  the  distance 
from  the  particle  to  the  target  is  increased,  such  that  the  an^e  subtended  by  the  target 
as  viewed  from  the  scattering  particle  becomes  small,  the  wavelength  dependence  also 
decreases.  This  result  can  be  anticipated  from  Eq.  (18)  as  dm  ax  0. 
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The  present  results  have  been  extended  to  the  case  of  the  attenuation  of  a laser 
beam  by  a natural  snowfall,  by  use  of  Eqs.  (18)  through  (33).  The  result  given  in  Fig. 
13  is  valid  for  both  a 0.6328-  or  a 10.6-jli  wavelength,  since  for  an  angle  of  0m  ax  = 
1.19  X 10‘3  deg  or  smaller,  the  difference  at  a liquid  water  content  of  0.25  g/m^  is 
less  than  0.82  percent. 

The  laser  testing  capability  developed  for  the  present  test  can  be  readily  expanded 
to  different  particulate  fields,  such  as  dust,  cloud,  or  rain.  The  technology  demonstrated 
for  a high  density  snowfall  of  limited  extent  may  also  be  applied  in  Range  G to  produce 
a longer  expanse  of  low-density  particulate  environment  for  laser  propagation  testing,  which 
more  nearly  duplicates  that  found  in  nature. 
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Figure  1.  Experimental  arrangement. 
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Figure  2.  Sketch  of  attenuation  field  geometry. 
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d/d  d/d 


Figure  3.  Typical  test  data. 
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Scattered 


Figure  4.  Scattering  by  single  particle. 
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Figure  5.  Radiation  intensity  distribution  on  detector. 
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Figure  7.  Signal  attenuation  ratio  for  single  particle. 
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Laser 

View 


Figure  8.  Typical  cell  for  analysis  of  environment. 
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Run  No . 9 


Figure  10.  Comparison  of 
resuits  • runs  9 


Run  No 


I 


10. 


and  predicted 


Experimental  Transmission 


Figure  11.  Comparison  of  calculated  and  experimental  transmission 
of  HeNe  laser  through  snowfield  - 0.9  sec. 
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Experimental  Transmission 


Figure  12.  Comparison  of  calculated  and  experimental  transmission 
of  CO2  laser  through  snowfield  - 0.9  sec. 
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Table  1.  Equipment  Description 


HeNe  Laser 


KeNe  Laser  Detector 


CO2  Laser 


CO2  Laser  Detector 


Spectra  Physics  Model  No.  124 

Power  Output  - 15  mw  continuous 

Beam  Diameter  - 1.5  mm 

Beam  Divergence  - 1 milliradian 

Metro! ogic  Radiometer  Model  No.  60-530 
Sensor  Head  - Silicon  photosensor  with 
radiometric  filter 
Sensor  Active  Area  - 0.385  cm2 
Frequency  Response  - 1 to  10  kHz 
Ranges  - 0.3  microwatts  to  30  milliwatts 

Sylvania  Corporation  941 P 

Power  Output  - 3 w continuous 
Beam  Diameter  - 4.5  mm 

Beam  Divergence  - 4 milliradians,  total  angle 

Molectron  Corporation  PE3-00 

Sensor  Element  - permanently  poled  pyroelectric 
crystal 

Sensor  Active  Area  - 1 x 1.2  mm 
Rise  time  - 5 psec 

Santa  Barbara  Research  Center  F256-2 
Sensor  Element  - Photoconductive 

Mercury  - Cadmium  - Telluride 
Sensor  Active  Area  - 1 x 1 mm 
Time  Constant  - <10  ysec 
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Table  2.  Test  Log 


Test 

No. 

5 

6 

7 

8 
9 

10 

n 

12 

13 

14 

15 

16 

17 

18 
19 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


Date  ' 
6/22/77 


6/28/77 


Snow 

Growth 

T1me> 

win 

8 

12 

5 

4 

3 

3 

4.5 

4 
4 

4 

5 

5 

4.5 

6 
6 

4 

3.2 

4.5 

5 

4.3 

3.5 
3.8 

3 

18 

6.5 

4 

5 


Dew 

Point 

66.8 


63.3 


36 

6 

60.8 

37 

6 

38 

6.5 

39 

6/28/77 

7 

60.8 

40 

7.5 

60.8 

41 

8 

62.7 

42 

10 

43 

9 

Remarks 

No  data  traces 
Noisy  traces 
Noise  on  HeNe  trace 


Zero  on  CO2  trace 
No  traces  recorded 

No  field  camera  picture 


No  pictures 
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Table  3.  Sample  Snowflake  Size  Distribution  Data 


Run  No.  = 19 

Time  = 

0.9  sec 

CELL  1 

CELL  2 

CELL  3 

Area 

Number 

Area 

Number 

Area 

Number 

(Picture  Points) 

of  Flakes 

(Picture  Points) 

of  Flakes 

(Picture  Points) 

of  Flakes 

0.11 500 

02 

2 

0.10860 

02 

7 

0.92000 

01 

5 

0.21 670 

02 

3 

0.29000 

02 

2 

0.25000 

02 

1 

0.29000 

02 

2 

0.30000 

02 

2 

0.30670 

02 

3 

0.40250 

02 

4 

0.47000 

02 

1 

0.47330 

02 

3 

0.62000 

02 

1 

0.60000 

02 

2 

0.0 

2 

0.73000 

02 

1 

0.69000 

02 

2 

0.75500 

02 

2 

0.94000 

02 

0 

0.84000 

02 

1 

0.91000 

02 

1 

0.26800 

03 

1 

0.0 

1 

0.24600 

03 

1 

0.23100 

03 

1 

0.18430 

03 

4 

0.29430 

03 

3 

0.37670 

03 

6 

0.44440 

03 

5 

0.39830 

03 

4 

0.66000 

03 

1 

0.64050 

03 

2 

0.57500 

03 

1 

0.89730 

03 

3 

0.79650 

03 

2 

0.94200 

03 

1 

0.13030 

04 

4 

0.14060 

04 

3 

0.12460 

04 

4 

0.17090 

04 

2 

0.19570 

04 

1 

0.20980 

04 

3 

0.26110 

04 

2 

0.27320 

04 

1 

0.0 

0 

0.40720 

04 

2 

0.0 

0 

0.38080 

04 

1 

0.0 

0 

0.0 

0 

0.0 

0 

0.0 

0 

0.74770 

04 

1 

0.0 

0 

210  picture  points  = 1 inch 


Table  4.  Attenuation  Data 


1/1 


Time-sec 


Time-sec 


Run 


No. 

0.74 

0.82 

0.9 

0.98 

1.06 

0.74 

0.82 

0.9 

0.98 

1.06 

7 

1.0 

0.75 

0.4 

0.5 

0.75 

- 

- 

- 

- 

- 

8 

1.0 

0.83 

0.43 

0.5 

0.78 

1.0 

0.86 

0.64 

0.62 

0.82 

9 

0.90 

0.88 

0.5 

0.53 

0.80 

1.0 

0.90 

0.64 

0.80 

0.96 

10 

0.80 

0.69 

0.40 

0.58 

0.69 

1.0 

0.98 

0.68 

0.80 

0.92 

12 

0.91 

0.69 

0.63 

0.77 

0.86 

1.0 

1.0 

0.70 

0.78 

0.92 

14 

1.0 

0.70 

0.58 

0.73 

0.85 

1.0 

0.84 

0.60 

0.80 

0.92 

15 

0.97 

0.92 

0.51 

0.74 

0.84 

1.0 

0.98 

0.44 

0.64 

0.86 

16 

0.98 

0.61 

0.49 

0.73 

0.93 

1.0 

0.84 

0.60 

0.82 

0.88 

17 

1.0 

0.60 

0.36 

0.53 

0.85 

1.0 

0.89 

0.58 

0.75 

0.87 

18 

0.81 

0.32 

0.30 

0.78 

CO 

1 

0.88 

0.46 

0.54 

0.74 

0.88 

19 

0.95 

0.56 

0.32 

0.56 

0.95 

1.0 

0.64 

0.56 

0.70 

0.82 
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Table  4.  Concluded 


Time-sec 


Run 

No. 

0.66 

0.74 

0.82 

0.90 

0.98 

1.06 

25 

1.0 

1.0 

0.95 

0.93 

0.84 

0.88 

26 

0.98 

0.91 

0.87 

0.71 

0.69 

0.76 

27 

1.0 

0.93 

0.89 

0.71 

0.64 

0.75 

28 

1.0 

1.0 

0.98 

0.83 

0.70 

0.80 

29 

1.0 

1.0 

1.0 

0.91 

0.84 

0.84 

30 

1.0 

0.91 

0.91 

0.87 

0.73 

0.73 

32 

0.98 

0.82 

0.23 

0.20 

0.31 

0.53 

33 

0.98 

0.96 

0.73 

0.51 

0.62 

0.77 

34 

0.97 

0.97 

0.92 

0.85 

0.71 

0.80 

35 

0.93 

0.91 

0.82 

0.65 

0.61 

0.75 

36 

0.94 

0.92 

0.88 

0.76 

0.63 

0.75 

37 

1.0 

1.0 

0.96 

0.75 

0.71 

0.82 

38 

0.81 

0.83 

0.82 

0.60 

0.60 

0.67 

39 

1.0 

1.0 

0.93 

0.63 

0.54 

0.65 

40 

0.80 

0.79 

0.68 

0.51 

0.53 

0.71 

41 

0.92 

0.92 

0.87 

0.55 

0.57 

0.77 

42 

1.0 

1.0 

0.58 

0.28 

0.54 

0.70 

43 

1.0 

1.0 

0.82 

0.45 

0.59 

0.70 

HeNe 


Time- sec 


1.14 

0.66 

0.74 

0.82 

0.90 

0.98 

1.06 

1.14 

0.95 

1.0 

1.0 

1.0 

0.98 

0.82 

0.82 

1.14 

0.91 

1.0 

1.0 

0.94 

0.87 

0.71 

0.82 

0.90 

1.0 

1.0 

1.0 

1.0 

0.81 

0.71 

0.81 

0.85 

1.0 

1.0 

1.0 

1.0 

0.95 

0.78 

0.88 

1.0 

0.91 

1.0 

1.0 

1.0 

0.95 

0.88 

0.88 

0.91 

0.80 

1.0 

1.0 

1.0 

1.0 

0.83 

0.79 

0.83 

0.71 

1.0 

0.91 

0.47 

0.20 

0.37 

0.63 

0.82 

0.83 

1.0 

1.0 

0.97 

0.70 

0.67 

0.82 

0.90 

0.84 

1.0 

1.0 

1.0 

1.0 

0.79 

0.89 

0.97 

0.90 

1.0 

1.0 

1.0 

0.82 

0.70 

0.82 

0.93 

0.83 

1.0 

1.0 

1.0 

0.90 

0.72 

0.81 

0.92 

0.92 

1.0 

1.0 

1.0 

0.82 

0.70 

0.82 

0.97 

0.72 

1.0 

1.0 

1.0 

0.80 

0.80 

0.77 

0.93 

0.80 

1.0 

1.0 

1.0 

0.82 

0.63 

0.74 

0.89 

0.-72 

1.0 

1.0 

1.0 

0.73 

0.68 

0.80 

0.91 

0.90 

1.0 

1.0 

1.0 

0.61 

0.60 

0.72 

0.87 

0.78 

1.0 

1.0 

0.78 

0.30 

0.56 

0.73 

0.85 

0.80 

1.0 

1.0 

1.0 

0.68 

0.60 

0.80 

0.91 
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COMPARISON  OF  EXPERIMENTAL 
AND  THEORETICAL  RESULTS 
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p/p, 


p/p, 


AED 


Time  f sec 


Run  No.  25 
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□ HeNe 
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_J I 
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Time , sec 


1.2 


Run  No . 26 
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Run  No.  33 


AEDC-' 


'd/d 


Run  No.  34 


Time,  sec 
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Run  No . 36 


Run  No.  37 
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p/p, 
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A 

Ad 

Ae 


Ae.i 

AIo.K 

Al 


Dd 


Di 


Dl 

Do 

Dp 

ADp.i 

D, 

I 

lo 

Is 

k 

Lc 

Lt 

i 


NOMENCLATURE 

Area 

Area  of  detector 

Effective  area  of  extinction  due  to  particles 

Effective  extinction  area  of  a single  i^’’  size  particle 

Area  which  receives  radiation  intensity  after  attenuation  by  K cells 

Total  extinction  area  of  the  cell 

Cross-sectional  area  of  laser  beam 

Effective  scattering  area  of  particle 

Detector  diameter 

Equivalent  melted  diameter,  mm 

Diameter  of  laser  beam  at  snowfield 

Equivalent  melted  value 

Particle  diameter 

Particle  size  range  increment 

Effective  scattering  diameter  = \/2  x Dp 
Intensity  of  radiation 

Intensity  of  incident  radiation  on  particle 

Scattered  radiation  intensity 

lirifK 

Dimension  of  cell  in  direction  along  laser  beam 

Total  extent  of  snowfield 

Distance  from  scattering  particle  to  detector 
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m liquid  water  content  of  cloud,  g/m^ 

Nj  Total  number  of  particles  in  the  i^^  size  range  per  unit  volume  per  unit 

size  range,  l/m^-mm 

Nl,!  Number  of  particles  in  the  i***  size  range  encountered  by  laser  beam  in 

a particular  cell 

^SV,i  Number  of  particles  in  side  view  of  ceU  in  the  i*^  size  range 

P Power  incident  on  detector  or  target 

Pq  Initial  power  in  laser  beam 

r Radial  coordinate  of  point  on  detector 

S Distance  into  cloud,  ft 

W Dimension  of  cell  used  for  snowfleld  calibration 

X IT  Dp/X 

a Constant  = 3.67 

d Scattering  angle 

dm  ax  Angular  extent  of  detector 

X Wavelength 

p Density  of  water  (taken  to  be  Ig/cm^) 

<j>  Angular  coordinate  of  point  on  detector 
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